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On the General Solution of Laplace’s Equation and the Equation 
of Wave Motions , and on an undulatory explanation of 
Gravity . By E. T. Whittaker, M.A. 

1. Laplace’s Equation. 

The first part of this communication is concerned with the 
well-known equation which was first introduced in Laplace’s 
memoir on the ring of Saturn , and which is now of such impor¬ 
tance in mathematical physics. The general solution has not 
hitherto been found in a finite form : we shall show that it is 
given in a fini te form by the following result : 

The general solution of Laplace’s equation 

9w , aw , aw 

-+-+ -= o 

a# 2 a2/ 2 a* 2 

,27T 

is Y =J f(x cos v + y sin v + iz, v) dv, 

where f is an arbitrary function of the two arguments x cos v + y 
sin v + iz and v. 

To establish this, let Y (x f y , z) be any solution (single-valued 
or many-valued) of Laplace’s equation. Since the differential 
equation is unaffected by a change of origin, and since any func¬ 
tion of the two arguments (x cos v-\-y sin v + iz) and v is also a 
function of the arguments (x' cos v + y ; sin v + iz r ) and v, where 
(1 x f , ?/, zf) are the coordinates of the point (x, y , z) referred to 
parallel axes through any other origin, it clearly will not affect 
the generality of the proof if any particular point be taken as 
origin ; we shall therefore suppose that the origin is taken at 
some point at which some branch of the function Y ( x , y , z) is 
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regular, so that this branch of Y ( x , y , 2) can be expanded in the 
neighbourhood of the origin in the form 

Y (x, y , 2) = + + + + + ^ + + 

+/ 2 ^+^ 3 + 

where a Q , a x , b lf ... . are constants, and the series is absolutely 
and uniformly convergent for a domain of finite extent. 

Substituting this expansion in Laplace’s equation, and equat¬ 
ing to zero the coefficients of the various powers of x , y, and z, 
we obtain an infinite number of linear relations between the 
constants a 0 , a T , b x , . . . There are \n(n— 1) of these relations 
between the %{n + i)(w + 2) coefficients of the terms of any 
degree n in the expansion of Y ; so that only {^(n + i){n+ 2) 
—^n(n—i)} or (2 n-\- 1) of the coefficients of terms of degree n in 
the expansion of Y are really independent. 

Moreover, it is clear that the terms of degree n in Y must 
satisfy Laplace’s equation, quite independently of the terms of 
other degrees ; and therefore the terms of degree n in Y must be 
a linear combination of (271+1) linearly independent solutions of 
degree n. 

To find these solutions, consider the expansion of the quantity 
(x cos v-\-y sin v + iz) n as a sum of sines and cosines of multiples 
of v, in the form 

(x cos v+y sin v + iz) n = g 0 {x , y, z)+g x (x, y, z) cos v 

+g 2 (x, y, 2) cos 2v+ .... +g n (x , y, z) cos nv 
+ h x (x, y , z) sin v + h 2 (xj y , z) sin 2V+ ... . 

+ h n (x, y , z) sin nv. 

The (271+ 1) quantities g (x, y , z) and h(x, y,z) are clearly 
linearly independent homogeneous polynomials of degree n in 
x, y, z ; and each of them must satisfy Laplace’s equation, since 
(x cos v-\-y sin v + iz) n does so. These are therefore the (2 n-\- 1) 
independent solutions of Laplace’s equation, of degree n. 

Now since by Fourier’s theorem we have the relations 

1 f 27r 

g m (Xj y, z) = - (x cos v-\-y sin v + iz) n cos mv dv 

1 f 27r 

h m (x, y, z) — - (x cos v + y sin v + iz) n sin mv dv, 

it follows that each of these (271+1) solutions can be expressed 
in the form 

* n2 ir 

f(x cos v+y sin v + iz, v) dv, 

and therefore any linear combination of these (2^+1) solutions 
can be expressed in this form. That is, the terms of degree n 
in Y ( x , y, z) can be expressed in this form, where n is any 
positive integer ; and therefore Y (x, y , z) (which was taken to 
be any solution of Laplace’s equation) is itself of this form ; 
which establishes the theorem. 
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2. The Partial Differential Equation of Wave-Motions. 

It can in the same way be shown that the general solution 
of the general differential equation of wave-motions, namely 

9 2 Y_j_ 9 2 Y_^ 9 2 Y_^ 2 9 2 Y 

9 sc 2 9 y 2 dz 2 dt 2 

as 

Y — j jy (as sin v cos ip + y sin v sin \p + z cos -u -}-&£, u, \p) dv dip, 

where f is an arbitrary function of the three arguments 

x sin v cos \p + y sin v sin \j/ + z cos v + ht, v, and ip. 

From this result it follows, by the application of Fourier’s 
double-integral theorem, that any solution of the equation 


9 2 Y + 9 2 Y + 9 2 Y 
dx 2 dy 2 dz 2 



can be analysed into simple plane waves. 


3. An TJndulatory Explanation of Gravity. 

In connection with the last property it can be shown that 
distributions of simple waves exist such that , although each indi¬ 
vidual wave is periodic with respect to the time, yet the total dis¬ 
turbance at any point does not vary with the time. The s im plest 
example of this is the following. Suppose a particle to be emit¬ 
ting spherical waves, such that the disturbance due to the waves 

whose wave-length lies between — and — 27r - is 

p p. + dp 

2 dp sin plfft—r) 

. , 

Tip r 

where Y is the velocity of propagation of the waves. Then 
after the waves have reached the point r, so that (Vt—r) 
is positive, the total disturbance there (due to the sum of all the 
waves) is 

^psin^Yi-r) , 
r-rJ 0 p 

or —f S * n ^ dy, where p(Vt — r)—y, 

TrrJ 0 y - 


i.e. the disturbance is independent of the time and is everywhere 
proportional to the gravitational potential of the particle. It is 
easily seen that for all such distributions the total disturbance 
can always be represented by a solution of Laplace 3 s equation. 

x x 2 
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Since Laplace’s equation is the equation satisfied by the com¬ 
ponents of any field of force due to gravitating bodies, this result 
suggests an explanation of the propagation of gravitation. For 
we thus see that it is possible to analyse the field of force due to 
a gravitating body, by a “ spectrum analysis ” as it were, into an 
infinite number of constituent fields; and although the whole 
field of force does not vary with the time, yet each of the con¬ 
stituent fields is of an undulatory character, i.e. the disturbance 
at any point reproduces itself after the lapse of any time at 
another point whose distance is proportional to the elapsed time. 
But this result assimilates the propagation of gravitation to that 
of light; for the undulatory phenomena just described, in which 
the varying vector is a gravitational force perpendicular to the 
wave-front, may be compared with those made familiar by the 
electromagnetic theory of light, in which the varying vectors 
consist of electric and magnetic forces parallel to the wave-front. 
The waves are in other respects of the same type in the two- 
cases, and it seems probable that an identical property of the 
medium ensures their transmission through space. This undula¬ 
tory theory of gravity would require that gravity should be 
propagated with a finite velocity, which, however, need not be the 
same as that of light, and may be enormously greater. 

Of course this investigation does not explain the cause of 
gravity; all that is done is to show that the propagation 
across space of forces which vary as the inverse square of the 
distance does not require for its explanation any other property of 
the medium than one which has long been known and accepted. 


On the Mean Distance of a Planet , as a Function of Three Helio¬ 
centric Distances and the Observed Times. By Shin Hirayama. 

§ i. Let r ZJ r 2 , r 3 be three heliocentric distances of a planet; 
t x , t 2 , t 3 the corresponding observed times ; 
a the mean distance. 

The radius vector of a planet must satisfy the differential- 
equation 

_ 1+2 =Q 

dt* Vs ® 

putting r 2 = s. 

If we let 

0 1 = k ( t 3 — t 2 ), d 2 = h { t 3 — t z ), 0 3 = hfa — tj ), log h = 8*235581 — 10,. 

and express the value r 2 in terms of t by expansion into series,, 
we have 

8 = r 2 = a Q -f a x t + a 2 t 2 + a 3 tf -f a 4 t\ 
where a Q a z , a 2 , a 3 , a 4 are five constants independent of time. 
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